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ABSTRACT 

We present initial results from a Chandra survey of a complete sample of the 8 nearest (z < 0.04) 
ultraluminous IR galaxies (ULIRGs), and also include the IR-luminous galaxy NGC 6240 for com- 
parison. In this paper we use the hard X-rays (2-8 kev) to search for the possible presence of an 
obscured AGN. In every case, a hard X-ray source is detected in the nuclear region. If we divide 
the sample according to the optical/IR spectroscopic classification (starburst vs. AGN), we find that 
the 5 "starburst" ULIRGs have hard X-ray luminosities about an order-of-magnitude smaller than 
the 3 "AGN" ULIRGs. NGC 6240 has an anomalously high hard X-ray luminosity compared to the 
"starburst" ULIRGs. The Fe Ka line is convincingly detected in only two ULIRGs. The weakness of 
the Fe-K emission in these ULIRGs generally suggests that the hard X-ray spectrum is not dominated 
by reflection from high Nh neutral material. The hard X-ray continuum flux ranges from a few x 10 -3 
to a few x 10~ 5 of the far-IR flux, similar to values in pure starbursts, and several orders-of-magnitudc 
smaller than in Compton-thin AGN. The upper limits on the ratio of the Fe Ka to far-IR flux are be- 
low the values measured in Compton-thick type 2 Seyfert galaxies. While very large column densities 
of molecular gas are observed in the nuclei of these galaxies, we find no evidence that the observed 
X-ray sources are obscured by Compton-thick material. Thus, our new hard X-ray data do not provide 
direct evidence that powerful "buried quasars" dominate the overall energetics of most ultraluminous 
infrared galaxies. 

Subject headings: evolution-galaxies:evolution-X-rays: galaxies-X-rays 



1. INTRODUCTION 

Ultraluminous IR galaxies (ULIRGs) are defined to 
be galaxies whose bolometric luminosity exceeds 1O 12 L 
and whose emis sion is dominated by the mid- and far- 
infrared bands ijSanders fe Mirabel 1 1996th 3 Thev ar e 
thought to be powered by either an AGN ( Sanders 1999) , 
a very powerful nuclear starburst l) Joseph! 1999|) . or both. 
It is possible that these galaxies are an evolutionary stage 
of normal galaxies in which a large amount of dense gas 
that has been driven to the nucleus is both feeding a nu- 
clear black hole and fueling a very high star-formation 
rate l|Norman fc ScovilleilT98cl . 

ULIRGs generically have highly disturbed morpholo- 
gies suggestive of severe tidal effects, and often have 
double n uclei suggestive of the late stages of a galaxy 
merg er l|Clements. Sutherland. McMaho n. fc Saunders I 
1996). If elliptical galaxies are the by-product of similar 
mergers, then ULI RGs should have be en very common in 
the early universe l)Oenzel et al.12 001'1. This is consistent 
with the properties of the sub-mm source population, 
and the corresponding estimates for the star-formation 
history of the universe l)Blain et al.lll999|) . 

ULIRGs offer ideal local laboratories for studying the 
processes occurring in the high-redshift sub-mm sources. 

3 We assume Hq = 70 km s -1 Mpc -1 . 



As such, decoupling the effects of an AGN and star- 
formation is important for assessing the role of ULIRG- 
like objects in building present-day elliptical galaxies, in 
accounting for the total amount of energy supplied by 
AGN over cosmic time, and in establishing the strong 
correlation between the mass of a supermassive black 
hole and the properties of its "host" galaxy s pheroid 
ijFerrarese fe MerrittJl2000l iTremaine et alJl2002|h 

The nature of the central power-source in ULIRGs re- 
mains uncertain, owing to its location behind large col- 
umn densities of obscuring material. This highlights the 
importance of observations in the hard X-ray, far-IR, and 
radio domain, which most effectively pierce the obscuring 
material. However, such observations to date have not 
provided decisive evidence to break the AGN/starburst 
deadlock (because they have lacked either the appropri- 
ate spatial resolution or have not provided unambigu- 
ous spectroscopic discriminants between AGN and star- 
bursts). 

The optical and near/mid-IR regimes con- 
tain many potentially important spectral di- 
agnostics. Such observations of ULIRGs vari- 
ously clas sify them a s Seyfer ts, LINERs or star- 
burst s iGenzel et alJ 119981 llmanishi fe Dudley ! 
20001 iKewlev. Heisler Dopita. fe Lumsdenl 120011 
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200C ; iVeilleux . Kim. Sanders. Mazzarella. fe Soifer I 

1995; IVeilleux. Kim, fc Sanded 

Veilleux. Sanders. fc Kimll999|) . If the ULIRGs optically 
classified as LINERs are considered to be starburst- 
driven, then there is good over all agreement between the 
optical and IR classifications jL utz. Veilleiix. fc Genzell 
119991 iTaniguchi. Yoshino. Ohvama. fc Nishiural 119991 , 
and in general nearby LINERs exhibiting only narrow 
Ha lines tend to have X-ray properties more consistent 
with a starburs t rath e r than an AGN interpretatio n 
l|Eracleous et al.l 12001 iTerashima. Ho. fe Ptakl l2000|) . 
However, X-ray observations of the infrared-luminous 
merger NGC 6240 (which appears to be starburst- 
dominated LINER from an optical/IR perspective) 
reveal a powerful AG N that may dominate its energetics 
l|Vignati et alJlT999^ . 

Observations of ULIRGs with the Chandra X-ray 
Observatory promise to be enlightening for this is- 
sue. Chandra has large effective area, high spatial 
resolution (FWHM ~ 1"), and is capable of moder- 
ate spectral resolution imaging spectroscopy. All of 
these may be necessary to determine the physical ori- 
gin of the hard X-ray emission previously detected from 
ULIRG s (iRisaliti. Gilli Maiolino fe Salvatil|200G{l Re- 
cently. iGallagher et alJ l)2002j) . iXia et all (|2002^ . and 
iClements et alJ l)2002|) have presented detailed analysis 
of Chandra observations of the ULIRGs Mkn 231, Mkn 
273, and Arp 220 respectively (see also iMcDowell et al.l 
2003). Our approach here is to examine a sample of 
ULIRGs that is large enough to document the properties 
of the class as-a-whole. 

The eight ULIRGs listed in Table 1 comprise both a 
flux-limited sample with Fqq^ > 10 Jy, and a volume- 
limited sample out to cz = 14000 km s" 1 (D < 200 Mpc), 
and NGC 6240 is included since it is an archetype of IR- 
luminous galaxie s shown by hard X-ray dat a to harbor 
a powerful AGN l)Iwasawa fc Comastrilll998l) . Note that 
the galaxies are all at a redshift of ~ 0.04 with the ex- 
ceptions of Arp 220 and NGC 6240 (z ~ 0.02). The joint 
analysis of XMM-Newton and Chandra data for two of 
these galaxies (IRAS 05189-2524 and UGC 05101) will be 
presented in future work. An analysis of the spatially- 
resolved soft X-ray emission and a comparison to models 
of galactic winds will also be presented elsewhere. 

The organization of the present paper is as follows. In 
section 2 we describe the data reduction and analysis 
techniques. In section 3 we discuss the results of spatial 
and spectral fitting of the data. In section 4 we discuss 
the implications of these results. 

2. DATA REDUCTION 

The observation dates and net exposure times for our 
targets are listed in Table 1. The data for Arp 220, 
Mkn 231, Mkn 273 and NGC 6240 were downloaded 
from the Chandra archive. The data reduction was per- 
formed using CIAO 2.2 and th e data were reprocesse d 
using CALDB 2.9. XAssist (see lPtak fc Griffiths! 12002]) . 
which is a software package that assists in data repro- 
cessing, initial source selection and analysis, was also 
used in the analysis of these data. Briefly, XAssist per- 
forms the basic data reduction steps recommended by 
the CXC "threads" (note that the 0.5 pixel position ran- 
domization is also removed, which is an optional step). 
Sources are detected using wavdetect, the background 



light curve is examined and time of background flaring 
are removed, and the image of each source is fitted with 
an elliptical gaussian model in order to determine the 
spatial extent (see section 3). The data were analyzed 
in full (0.3-8.0 keV), soft (0.5-2.0 keV) and hard (2.0- 
8.0 keV) bandpasses. For comparison with previous X- 
ray studies a region was selected (manually) that encom- 
passed the entire X-ray extent of the galaxy (determined 
from the full-band images) and spectra were extracted 
for these "global" regions. The spectra were binned to 
20 counts/channel to allow use of the x 2 statistic. 

Figure 1 shows adaptively-smoothed soft and hard 
band images, with the global and nuclear regions marked. 
We show the soft band images here for completeness, but 
defer their analysis. Note that the nuclear regions were 
chosen to maximize the signal-to-noise of the nuclear 
spectrum while the emphasis for the global regions was 
complete coverage of the X-ray emission of the galaxies 
(in part for comparison with previous X-ray studies). We 
checked for counterparts to all X-ray sources within the 
central 8' of the observation (where the PSF size is < 1") 
using HEASARC databases containing stars (USNO and 
GSC2.2), QSOs, 2Mass sources and the FIRST survey. 
In general, only a few counterparts in each field were 
found. These show that the astrometric solutions are 
good to ~ 0.5 — 1", but could not be used to apply reli- 
able astrometric corrections. 

3. SPATIAL ANALYSIS 

3.1. Nuclear Source Sizes 

The extents derived from the spatial fitting procedure 
are listed in Tabl e 2. The fi tting was performed using 
the C fit statistic (|Cashlll9'79Tl and the errors correspond 
to AC = 4.605 (90% confidence interval for 2 interesting 
parameters). In all but three cases, the nuclear 2.0-8.0 
keV source is consistent with being unresolved (FWHM 
< 1"). The hard X-ray emission in the nuclei of both 
IRAS 17208-0014 and Arp 220 are resolved and asym- 
metric, with a major axis FWHM of 2.6 " (0.9 kpc) in 
Arp 220 and 4.7 " (3.9 kpc) in IRAS 17208-0014. Point 
sources are present in the FOV with Gaussian fit FWHM 
parameters of ~ 1", indicating that the elongation in 
IRAS 17208-0014 is not due to errors in the aspect solu- 
tion. In NGC 6240, the hard X-ray nuclear emission is 
resolved into two sources coincident with the radio nuclei, 
and diffuse emissio n, consistent with the full-ban d HRC 
results reported in iLira. Ward. Zezas fc Murray! ( 2002) 
and the ACIS results reported in lKomossa et alJ f2003). 
However for consistency with the rest of the sample we 
model the hard X-ray flux with a single Gaussian to de- 
termine the extent of the total (i.e., combined) nuclear 
flux, noting that the individual nuclei would not have 
been resolved at redshifts > 0.04 (see also below). 

In order to further investigate the presence of unre- 
solved hard X-ray emission the 2-8 keV nuclear emission 
was fit with a Gaussian model convolved with a PSF im- 
age generated using the CIAO tool mkpsf (using a pho- 
ton energy of 3 keV, where the 2-8 keV photon spectrum 
generally peaks for these sources). For the galaxies with 
a nuclear point source, the main impact of including the 
PSF was to reduce the source extent parameter a from 
~ 0.5"to ~ 0.25": that is, the model PSF convolved with 
a Gaussian with er=0.25"(FWHM = 0.6 arcsec) provides 
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Table 1. Basic data for the ULIRG sample. 



Galaxy 


Position 


z 


ocale 


Date 


Exposure 


galactic 




(J 2000) 




(kpc/ ) 




(ks) 


/-iri20 — 2*\ 
(10 cm J 


IRAS05189-2524 


05 21 01.5 -25 21 45 


0.0426 


0.82 


10/30/2001 
01/30/2002 


6.6 
14.8 


4 


UGC 05101 


09 35 51.6 +61 21 11 


0.0394 


0.76 


05/28/2001 


49.3 


3 


Mkn 231 


12 56 14.2 +56 52 25 


0.0422 


0.82 


10/19/2000 


34.5 


1 


Mkn 273 


13 44 42.1 +55 53 13 


0.0378 


0.74 


04/19/2000 


44.2 


1 


Arp 220 


15 34 57.1 +23 30 11 


0.0181 


0.36 


06/24/2000 


56.4 


4 


NGC 6240 


16 52 58.9 +02 24 03 


0.0245 


0.49 


07/29/2001 


36.6 


6 


IRAS17208-0014 


17 23 21.9 -00 17 00 


0.0428 


0.83 


10/24/2001 


48.5 


10 


IRAS20551-4250 


20 58 26.9 -42 39 


0.0428 


0.83 


10/31/2001 


40.4 


4 


IRAS23128-5919 


23 15 47.0 -59 03 17 


0.0446 


0.86 


09/30/2001 


33.8 


3 



Note. — Positions and rcdshifts were obtained from NED. Galactic Nh values are determined from Dickey & Lockman (1990) using the 
LHEASOFT tool "nh" (http://hcasarc.gsfc.nasa.govi 



Table 2. Nuclear Spatial Properties 



Galaxy 




Soft Emission 








Hard Emission 






(") 


°v (") 





0~x (' 


") 


°v (") 


e 


Arp 220 


1.05 (1.03-1.06) 


0.97 (0.95-0.98) 


64 (59-67) 


0.99 


(0.86-1.14) 


0.68 (0.57-0.82) 


190 (175 - 210) 


IRAS05189-2524 


0.46 (0.42-0.49) 


=C T x 




0.44 


(0.41-0.47) 


= 0~x 




IRAS05189-2524* 


0.50 (0.46-0.53) 


=C T x 




0.46 


(0.44-0.48) 


=<J X 




IRAS17208-0014 


2.02 (1.62-2.53) 


1.45 (1.20-1.76) 


57 (34-81) 


1.05 


(0.89-1.26) 


0.58 (0.48-0.74) 


66 (54-81) 


IRAS20551-4250 


1.53 (1.15-1.86) 


1.06 (0.76-1.33) 


150 (127-170) 


0.65 


(0.53-0.84) 


=o~x 


IRAS23128-5919 


1.35 (0.94-1.70) 


0.77 (0.55-1.06) 


146 (131-168) 


0.43 


(0.38-0.48) 


=<?x 




Mkn 231 


0.43 (0.40-46) 


=C T x 




0.42 


(0.40-0.44) 


= 0~x 




Mkn 273 


1.17 (1.06-1.29) 


0.61 (0.53-0.69) 


45 (38-52) 


0.59 


(0.55-0.64) 


0.47 (0.43-0.50) 


141 (128-152) 


NGC 6240 


1.85 (1.72-1.99) 


2.79 (2.58-3.04) 


158 (150-166) 


0.80 


(0.72-0.88) 


1.19 (1.12-1.26) 


14.5 (8.9-21.4) 


UGC 05101 


0.97 (0.83-1.10) 


=&x 




0.55 


(0.46-0.63) 


= °~x 





Note. — Fit results based on fitting an elliptical gaussian model where surface brightness — Ne ' 2cTx 2<T y , and x, y arc rotated from the 
right ascension axis by 9 (given in degrees). Note that these fits were performed only to the nuclear source in images extracted from the central 
5"of these galaxies, and as shown in Figure 1 often emission is present on larger spatial scales. 
*2nd observation of IRAS 05189-2524. 



a good fit to the nuclear point sources in these ULIRGs. 
We regard this as an upper limit to the true source size, 
since it may arise from errors in the aspect solution. 

As stated above, NGC 6240 is resolved into two faint 
sources in addition to diffuse flux. Modeling this flux as 
a combination of two point sources in addition to two 
elliptical Gaussians (to account for diffuse flux associate 
with each nuclei) we find that the northern and south- 
ern nuclei contribute 17% and 28% of the 2-8 keV counts 
from the central 9" (~ 0.5 kpc) of NGC 6240. In the 
cases of Arp 220 and IRAS 17208-0014, the diffuse emis- 
sion was modeled with elliptical Gaussians and the PSF 
x Gaussian model with a = 0.25" was placed at the cen- 
troid of the diffuse flux. This resulted in upper limits 
(for AC = 4.605) of 1.2 x 10~ 3 and 4.3 x 10~ 4 counts 
s~ 4 for the 2-8 keV count rate of a point source, respec- 
tively, or 2-10 keV fluxes of 4.5 x 10~ 14 erg cm -2 s _1 and 
8.3 x 10~ 15 erg cm" 2 s~ 4 (L 2 -io keV ~ 4 x 10 40 erg s" 1 
in both cases). Finally, we note that in Arp 220 there is 
a ~ 1.7" (0.6 kpc) offset between the hard and soft nu- 
clear positions, with the hard nucleus bei ng coincident 
with the western 6 cm nucleus discussed in lNorrisl (1988) 
and the soft nucleus being offs et to the NW of the hard 
nucleus ijClements et al.ll2002|) . The soft nucleus is not 
coincident with either radio nucleus and is most likely 



due to hot ISM observed through a region with relatively 
low extinction. Note that the optical nucleus of Arp 220 
is similarly offset by ~ 2" to the NW of the western ra- 
dio nucleus, which has likewise been interpreted as the 
result of a "hole" in the extinction toward the nucleus 
ijArribas. Colina fc Clements!l2001(l . 

3.2. Extra-nuclear Point Sources 

Table [3] lists the positions and luminosities of unre- 
solved sources that were detected within the X-ray ex- 
tent of the galaxies. The luminosities were computed 
assuming a power-law model with a photon index of 1.8. 
If these sources are associated with the corresponding 
host galaxy, then they have luminosities well in excess 
of the Eddington luminosity of a solar-mass black hole 
or neutron star X-ray binary, and hence are "ultralumi- 
nous X-ray objects" (ULXs; also known as intermediate- 
luminosity X-ray objects or IXOs). Note that the ULX 
in UGC 05101 has a very hard spectrum (i.e., it is 
not present in the soft-band image in Figure 1) which 
may be indicative of a highly-absorbed spectrum with 
Nh > 10 22 ~ 23 cm~ 2 . Since this source is not in the nu- 
cleus of UGC 05101, this would imply a large amount 
of int ernal absorption for this source. iGallagher et alJ 
(2002) also discuss an ULX in Mkn 231 which lies be- 
low our detection threshold of ~ 10 counts. Finally, the 



Arp 220 0.5-2.0 keV 



Arp 220 2.0-8.0 keV 




IRAS 051 89-2524 2nd Observation, 2.0-8.0 keV 



IRAS 1 7208-001 4 0.5-2.0 keV IRAS 17208-0014 2.0-8.0 keV 




Fig. 1 — 



ULX in IRAS 20511-4250 has an estimated hard X-ray 
luminosity that is ~ 20% of the luminosity of the nuclear 
source in this galaxy. However, this luminosity is highly 
uncertain given the low count rate of the source. The 2-8 
keV count rate of the ULX is - 25% of global 2-8 keV 
count rate. 




IRAS 20551-4250 0.5-2.0 keV 



IRAS 20551-4250 2.0-8.0 keV 




Fig. 1 



— (cont.) 
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Fig. 1. — (cont.) 0.5-2.0 keV and 2.0-8.0 keV images of the ULIRG sample. The outer and inner ellipses show the global and nuclear 
source region, respectively. The images were adaptively smoothed to a significance level of 2.8<x and scaled logarithmically. N is up and E 
is to the left in these images. 
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3.3. Large-Scale Emission 

In addition to the nuclear source and the extranuclear 
point sources, extended hard X-ray emission is often de- 
tected beyond the nuclear region on scales of ~ 10 kpc. 
The fraction of total hard X-ray flux in this spatially- 
extended component is < 10% in four cases (Arp 220, 
IRAS 05189-2524, Mkn 273, and Mkn 231), but ranges 
from 30 to 50% in the other cases (Figure 1). Note that 
in the case of IRAS 20551-4250, the extra-nuclear unre- 
solved source may account for up to ~ 50% of the extra- 
nuclear hard X-ray flux. 

4. SPECTRAL ANALYSIS 

As stated above, in every galaxy a nuclear core was 
detected in the 2-8 keV data. We proceeded initially by 
using XSPEC to fit these spectra with a simple absorbed 
power-law model. The results of this procedure are given 
in Table 4. Errors (corresponding to A% 2 = 4.605) are 
only derived for fits in which x 2 /dof < 1.5. Figure 2 
shows the data and model for these fits. These fits result 
in x 2 /dof < 1.5 only in the cases of IRAS 17208-0014 
and Mkn 231. 

We ne xt tried to fit each spectrum wi t h a " mekaF' 
plasma l|Liedahl. Osterheld. fc Goldsteinl ^995) plus 
power-law model. In these fits, absorption is applied 
to the entire spectrum with a hard lower-limit on the 
column density set to be 50% of the Galactic column, 
and a second absorber at the redshift of the ULIRG is 
applied only to the power-law component. Based on X- 
ray spectroscopy of starburst galaxies, the temperature 
of the plasma was limited to the range 0.3-2.0 keV. If 
the temperature was not constrained by the data then 
it was held fixed at 0.7 keV. Likewise, the abundance 
was limited to the range of Z = — 5Zq, and fixed at 
I.OZq if not constrained. The results are given in Ta- 
ble 5 and plotted in Figure 3. Here all fits have x 2 /dof 
< 1.5 except in the cases of Mkn 273, UGC 05101 and 
NGC 6240 (due to Fe-K emission, see below). Note that 
this more complex model is statistically preferred to the 
simple power- law model for Mkn 231. Significant residu- 
als are evident in the case of Arp 220 around 0.8 and 1.8 
keV, possibly the result of non-solar O and Si abundances 
or additional thermal components at different tempera- 
tures being required. A degeneracy between these sce- 
narios is often observed in nearby sta rburst galaxies (e.g., 
iWeaver. Heckman. fc Dahlemll2000|) . 

Chandra obser vations of br ight point sources may suf- 
fer from pile-up l)Davis) ;2001'l. This is only a concern for 
IRAS 05189-2524, the brightest ULIRG in our sample 
and the only one that is completely unresolved at both 
soft and hard X-ray energies (see also the pile-up d iscus- 
sion concerning Mkn 231 in lGallagher et al.ll2"002j) . We 
assessed the amount of pile-up by modeling the 2nd ob- 
servation, which had the higher count rate, using MARX. 
This analysis indicated that ~ 20% of the Chandra 
counts are piled-up in that observation. We checked the 
impact by fitting the two spectra from the two observa- 
tions of IRAS 05189-2524 with the plasma plus power-law 
model, additionally includin g the "pile-u p" convolution 
model in XSPEC (based on lDavisfl200l . We allowed 
the overall normalization and the pile-up grade migra- 
tion parameter to vary independently. This resulted in 
a power-law slope of T = 1.1, consistent with the indi- 



vidual fits, and showed that we are underestimating the 
2-10 keV flux of IRAS 05189-2524 by - 15%. Therefore 
pile-up does not affect our conclusions. 

4.1. Fe-Ka Emission 

In Mkn 273 and NGC 6240, strong Fe-K emission is 
present, and Fe-K emission is marginally detected in 
UGC 05101 (joint analysis of the Chandra and XMM- 
Newton data, to be presented elsewhere, verifies the ex- 
istence of the line) . In Table 6 we show the results of in- 
cluding a Gaussian line in the plasma + power-law model 
fits for IRAS 05189-2524, Mkn 231, Mkn 273, NGC 6240, 
and UGC 05101 (insufficient counts are detected from 
the other ULIRGs for these fits to result in useful con- 
straints). We also calculate the upper limits for a narrow 
line (physical width = 0.01 keV; the resolution of ACIS-S 
at ~ 6 keV is ~ 0.1 keV), and lines with physical widths 
of 0.1 and 0.5 keV at 6.7 keV using absorbed power-law 
models applied to the 2-8 keV spectra of IRAS 05189- 
2524 and Mkn 231. 6.7 keV would correspond to the Ka 
emission from Fe XXV. Such a line would be expected 
from an ionized scattering region, and we computed the 
upper limits assuming a non-negligible physical width 
to demonstrate the impact of a range of ionization states 
being present. N ote that the Fe-K line fr om IRAS 05189- 
2524 reported in lSevergnini et all l)2001|) is not detected 
in our Chandra (or XMM-Newton) data, although the 
significance of the line in the prior observations was low 
and our upper-limits are consistent with the errors. 

The hard X-ray continuum is usually very weak, and 
in many cases the resulting uncertainties on the contin- 
uum fit are substantial. Thus, a more direct (model- 
independent) upper limit to the Fe-Ka flux was deter- 
mined by interpolating the local background plus con- 
tinuum flux across the region of the Fe-Ka line in the 
spectra, and then setting an upper bound to the excess 
of line photons. As a consistency check, we computed 
the flux of the Fe-K line in Mkn 273 directly from the 
spectrum in this fashion, which resulted in a value of 
4.5 x 10 -14 erg cm" 2 s" 1 which is consistent with the 
value of 6.0 x 10 -14 erg cm -2 s _1 obtained from the fit 
listed in Table 6. The resulting upper limits on the Ka 
fluxes are listed in Table 7, where line fluxes derived from 
spectral fits are given for the galaxies listed in Table 6 
(i.e., those for which there were sufficient counts for spec- 
tral fitting in the Fe-K region). Note that the lines are 
typically about two orders-of-magnitude smaller than the 
flux of the strong Ka line in NGC 6240. 

Because the hard X-ray continua are so weak, the up- 
per limits on the equivalent widths (EW) of the Ka lines 
implied by the above are usually of-order a keV. The 
spectra of IRAS 05189-2524 and Mkn 231 have equiva- 
lent width upper limits of 0.1 keV and 0.5 respectively 
for a narrow Ka line (a = 0.01 keV) at 6.4 keV. Note 
that in the cases of UGC 05101, IRAS 17208-0014 and 
IRAS 20551-4250, insufficient counts exist in the 3-8 keV 
bandpass to even fit the power-law plus Gaussian model. 

In general, obscured (narrow-line) AGN in which 
the hard X-ray emission is dominated by reflec- 
tion/reproces sing tend to have l ines with EWs in excess 
of 1 keV fe.g.. lPtak et al.l l)1996l) and references therein), 
and similarly hot gas is expected to have Fe-Ka emis- 
sion at 6.7 keV with an equivalent width of ~ 0.6(Z/Zq) 
keV l|Rothernflug fc Arnaudlll985l) . Thus, in most cases 
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Table 3. Extra-nuclear. Point Sources 



Galaxy 


Position 


Offset 


Counts 


F2-10 keV 


L2-W kcV 


Notes 






("/kpc) 




10~ 15 erg cm~ 2 s _1 


10 40 erg s" 1 




Arp 220 


15 34 57.0 23 30 5.6 


7.3/2.8 


29.2 


1.8 


0.13 




IRAS 20551-4250 


20 58 26.5 -42 39 8.6 


8.4/7.5 


87.6 


6.9 


2.8 




IRAS 23128-5919 


23 15 46.7 -59 3 11.0 


4.7/4.0 


75.5 


6.4 


2.9 


Northern nucleus 


UGC 05101 


9 35 52.7 61 21 12.3 


8.3/6.8 


17.5 


1.2 


0.42 






Arp 220 



2 

E (keV) 



IRAS 05189-2524 - 




E (keV) 



Fig. 2.— 



RAS 05189-2524 (2nd Obs.) - 




E (keV) 



IRAS 17208-0014 



2 

E (keV) 



the upper bounds on the EW are not terribly restrictive. The upper limits on the Ka fluxes are more constraining 
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in these cases (as we will discuss below). 

4.2. More Complex Models 

The reduced \ 2 values of plasma + power-law + Fe- 
K line fits to Mkn 273, NGC 6240, and UGC 05101 
are > 1.5. In the case of UGC 05101 the fit is nev- 
ertheless statistically acceptable due to the low number 
of degrees of freedom. We added an addition plasma 
model (mekal) model to the Mkn 273 and NGC 6240 
fits. In the case of Mkn 273 this results in an accept- 



able fit with x 2 /dof = 41.2/40, however the photon in- 
dex of the power-law component is poorly constrained 
r = 1.7 (0.3 — 3.8). Since a motivation for these spectral 
fits is to compute accurate fluxes for comparison with 
global (i.e., spatially-averaged) spectral fits (see below), 
we fixed the photon index at 1.8, and the results of this 
fit are shown in Table 0D For NGC 6240 a more com- 
plex model was necessary in order to fit both the nu- 
clear and global spectru m. This model, mot ivated by 
the BeppoSAX fit given in lVignati et alJ l|1999|) . contains 
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Table 4. Power-law Fits to Nuclear Spectra 



Galaxy N H (10 22 cm' 2 ) T X 2 /dof 



Arp 220 0.0 0.9 39.5/21 

IRAS 05189-2524 0.0 -0.4 127/28 

IRAS 05189-2524* 0.0 -0.7 430.2/78 

IRAS 17208-0014 0.3 (< 0.6) 1.7 (1.0-2.3) 10.5/9 

IRAS 20551-4250 0.1 2.0 30.9/12 

IRAS 23128-5919 0.0 0.4 20.4/11 

Mkn 231 0.01 (<0.02) 0.44 (0.33-0.53) 57.0/53 

Mkn 273 0.0 0.3 302/49 

UGC 05101 0.0 1.0 35.5/12 

NGC 6240 0.39 1.4 431/136 



2nd observation 



two plasma components, a power-law (representing AGN 
emission scattered by a scattering region), an Fe-K com- 
plex (modeled with multiple Ga ussian components) and 
a Com pton reflection component ijMagdziarz fc Zgziarskil 
1995). See Table|5]for the fit parameters. We also found 
it necessary to include an additional Gaussian component 
with an observed (rest-frame) energy of ~ 2.2 (~ 2.3) 
keV, however note that the observed energy is close to 
the energies of significant absorption edges in the Chan- 
dra mirror response where the calibration tends to be the 
most uncertain, and accordingly this feature may not be 
real. 

4.3. Spatially-Averaged Spectra 

In order to determine the amount of flux extended be- 
yond the nuclear region we fitted the spectra extracted 
from the "global" regions shown in Figure 1 with a mod- 
els consisting of a plasma plus absorbed power-law (i.e., 



as fit to the nuclear spectra) and two plasmas and an ab- 
sorbed power-law. The parameters from the statistically- 
preferred model (based on the f-test statistic) are given 
in Table ITU1 The spectra from Mkn 273 and NGC 6240 
were fit with the same models used for the nuclear spec- 
tra, with the results shown in Tables |S1 and |5| 

4.4. Continuum Fluxes and Luminosities 

The observed fluxes and luminosities in the 0.5-2.0 keV 
and 2.0-10.0 keV bandpasses are given in Table ITTl de- 
rived from the plasma + power-law fits to the nuclear 
spectra, except for Mkn 273, UGC 05101 and NGC 6240 
where the fit in Table 6 was used. We also include the 
fluxes derived from fits to the "global" spectra (using 
the models shown in Tables 110(1 . from which it can be 
seen that on average ~ 50% and 75% of the 0.5-2.0 keV 
and 2.0-10.0 keV total flux originates in the ~ 0.5 — 1.0 
kpc scale nuclear regions. Since the sources with higher 
signal-to-noise require more complex models to fit their 
spectra it is likely that most or all ULIRGs are simi- 
larly complex (se e, e.g., XMM-Newton spectral fits in 
iBraito et al.ll2002T) . Fitting the spectra with models that 
are not sufficiently complex physically but are acceptable 
statistically may result in incorrect flux estimates, espe- 
cially in the case of the 2-10 keV flux estimates when 
the power-law slope is poorly constrained. In order to 
derive a model-independent estimate of the amount of 
extra- nuclear flux we also list in Table ITTl the ratio of the 
nuclear and global count rates in the 0.5-2.0 and 2.0-8.0 
keV bandpasses (with the 2.0-8.0 keV count rate ratios 
being listed with the 2-10 keV flux ratios). As expected 
the flux and count rates ratio are in fairly good agree- 
ment, but differ by up to ~ 20%. 

5. VARIABILITY 
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Fig. 3. 



Significant short-term variability in the 2.0-8.0 keV 
bandpass is not observed in the events extracted from the 
nuclear regions when tested by comparing the source and 
background photon arrival times with the K-S statistic , 
with the exception of Mkn 231 ijGallaeher et alJl2002j) . 
In Figure 4 we plot the long-term 2-10 keV light curves 
of the ULIRGs for which ASCA and/or BeppoSAX data 
exist (listed in Table I12f> . The Chandra fluxes where 
derived from fits to the "global" regions shown in Fig- 
ure 1, while the typical ASCA and BeppoSAX source 



regions were 4-6'. For the ASCA observation of UGC 
05101 and the two BeppoSAX observation of Arp 220, 
the source was not detected with sufficient counts for 
spectral modeling so the count rate was determined from 
the images and 2-10 keV fluxes were derived assuming the 
best-fitting spectral model to the Chandra data. While 
there is some apparent variability in IRAS 17208-0014 
and IRAS 20551-4250, we note that the total flux of all 
sources detected within the central 4' of the Chandra 
observations of these galaxies was 1.6xl0 -13 erg cm~ 2 
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Fig. 3.— (cont.) 



s _1 and 4.6xl0 -13 erg cm -2 s -1 , respectively. This 
suggests that some or all of the observed variability is 
due to source confusion in the ASCA data. In summary, 
hard X-ray variability is clearly observed in IRAS 05189- 
2524, Mkn 231, and Mkn 273 (the AGN ULIRGs) and 
in NGC 6240 at levels of ~ 40% or greater, while in the 
remaining ULIRGs no variability is observed that is not 
possibly due to aperture effects. 

6. DISCUSSION 



Our overall goal in this paper is to better understand 
the role of AGN in the ULIRG phenomenon. We can cast 
this in terms of three related questions: 1) Is the signa- 
ture of an AGN present in the Chandra hard X-ray data? 
2) Does an AGN dominate the hard X-ray emission? 3) 
Do the X-ray data imply that an obscured AGN dom- 
inates the bolometric energy output? To address these 
questions, we will consider the structure/morphology of 
the X-ray emission (information that Chandra uniquely 
provides), the spectral properties of the emission, and the 
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UGC 05101 



2 



E (keV) 




5 
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Fig. 3. — Plasma + power-law fits to the nuclear spectra of the ULIRG galaxies. The top panel shows the data and model and the 
data/model ratio is shown in the bottom panel. The plasma component model is shown with a dashed red line and the power-law component 
model is shown with a dot-dashed blue line. 



Table 5. Plasma + Power- law Fits to Nuclear Spectra 



Galaxy 


N H 

(10 22 cm- 2 ) 


kT (keV) 
(keV) 


Z/Z Q 


K** 


N H ,2 

(10 22 cm' 2 ) 


r 


X 2 /dof 


Arp 220 


0.6 (< 1.1) 


0.8 (0.6-1.1) 


2.1 (>0.2) 


1.3 (0.1 - 16) 


0.3 (<2.5) 


1.1 (0.3-1.8) 


19.5/17 


IRAS 05189-2524 


0.01 (< 0.3) 


0.7 (<1.3) 


0.01 (< 0.05) 


24. (15.-250) 


3.8 (2.2-6.0) 


1.0 (0.4-1.7) 


22.0/24 


IRAS 05189-2524* 


0.06 (< 0.2) 


0.4 (<0.8) 


0.00 (< 0.03) 


43. (16.-140) 


4.5 (3.4-5.8) 


1.1 (0.7-1.5) 


87/74 


IRAS 17208-0014 


1.0 (< 1.5) 


0.7(f) 


1.0(f) 


4.5 (<13.) 


0.6 (<39) 


1.9 (> 0.5) 


6.8/7 


IRAS 20551-4250 


0.6 (<0.8) 


0.3 (< 1.0) 


1.0(f) 


8.2 (0.6-21.) 


0.0 (<4.0) 


1.6 (0.8-3.8) 


11.2/9 


IRAS 23128-5919 


0.2 (< 0.5) 


0.7(f) 


0.0 (< 0.1) 


9.0 (2.5-18) 


4.8 (0.9-13.2) 


2.2 (0.5-4.6) 


8.0/8 


Mkn 231 


0.0 < 0.6) 


0.7(f) 


0.15 (0.02 - 3.2) 


1.8 (0.2-16) 


0.0 (< 0.6) 


0.2 (0.0-0.5) 


41.6/50 


Mkn 273 


0.3 


1.3 


0.06 


17 


24 


0.8 


72.6/45 


UGC 05101 


0.2 


0.4 


1.0(f) 


5.6 


0.3 


1.2 


29.3/9 


NGC 6240 


1.15 


0.65 


0.58 


300. 


0. 


0.80 


321/132 



Note. — Errors were determined using Ax 2 — 4.605 for fits with \ 2 jdof < 1.5. Parameters marked with "f" were frozen at the value shown. 
*2nd observation 
**Plasma model normalization in units of 10~ 9 - n ^"J^.> dV , n e = 

electron density in cm -3 , n p = Hydrogen density in cm -3 , D = 
luminosity distance to source in cm. 



ratio of the hard X-ray to bolometric (IR) luminosity. 

As noted in section 1, optical and IR spectroscopy of 
ULIRGs reach a surprising degree of agreement as to 
their classification. ULIRGs that are optically classified 
as Hll-regions or LINERs (type 1 or type 2 Seyferts) are 
nearly always classified a s starbursts (AGN) on the basis 
of mid-IR spectroscopy (iLutz. Veilleux. & Gen"zeil [i999: 
Tani guchi. Yoshino. Ohvama. fc Nishiural 119991) . Near- 
IR spectroscopy is also available for much of the sam- 
ple, and provides further diagnostic information, such 
as identifying AG N on the basis of detection of broad 
Paa or [Si VI] ijVeilleux. Sanders. Kr. Kiml IT999T) and 
the strength of PAH emi ssion and absorption features 
Ijlmanishi fc Dudley! l2000|) . Finally, AGN can be dis- 



criminated from starbursts (on an empirical, statisti- 
cal basis) using the F2B;.m/ftn ip, flux ratio, which is 
a me asure of dust temperature ijde Griip. Milev. fc Lubl 
I1987|) . This information is summarized in Table fTTH On 
the basis of these diagnostics, 5 ULIRGs are classified 
as starbursts (Arp 220, UGC 5101, IRAS 17208-0014, 
IRAS 20551-4250, and IRAS 23128-5919), and the other 
three are classified as AGN (Mkn 231, Mkn 273, and 
IRAS 05189- 2524). The analysis of a 3-4qm sp ectrum 
presented bv llmanishi. Dudley fc Malonevl l)2001|) shows 
that a highly-obscured AGN is likely to be present in 
UGC 05101. In what follows below, we will refer to 
these two classes as "starburst-ULIRGs" and "AGN- 
ULIRGs" respectively (for consistency we place UGC 
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Table 6. Fe-K Fits to Nuclear. Spectra 



v i < x, i aA v 


Bandpass 


Nu 


r 


Line E 


Line a 


Line EW 


v 2 /dof 




fkeVI 
l Kev J 


1022 cm — 2 




I'kpV'l 
l Kev i 


V Kev ) 






IRAS 05189* 


0.4-8.0 


4.5 (3.4-5.8) 


1.1 (0.7-1.5) 


6.4f 


O.Olf 


0.0 (<0.1) 


87/73 




2.0-8.0 


4.4 (2.4-6.7) 


1.1 (0.6-1.7) 


6.4f 


O.Olf 


0.0 (<0.1) 


73/57 




2.0-8.0 


4.8 (2.7-7.2) 


1.3 (0.7-1.9) 


6.7f 


O.lf 


0.1 (<0.4) 


71/57 




2.0-8.0 


4.7 (2.5-7.6) 


1.2 (0.6-2.1) 


6.7f 


0.5f 


0.1 (<0.7) 


73/57 


Mkn 231 


0.4-8.0 


(<0.8) 


0.2 (0.0-0.7) 


6.4f 


O.Olf 


0.2 (<0.5) 


38/49 




2.0-8.0 


0.0 (<1.6) 


0.3 (0.0-0.9) 


6.4f 


O.Olf 


0.2 (<0.6) 


12/24 




2.0-8.0 


0.0 (<1.5) 


0.3 (-0.1-0.9) 


6.7f 


O.lf 


0.1 (<0.5) 


16/24 




2.0-8.0 


0.0 (<2.3) 


0.5 (0.0-1.3) 


6.7f 


0.5f 


0.6 (<1.5) 


13/24 


Mkn 273 


0.4-8.0 


24 (9-40) 


0.9 (-0.5-2.6) 


6.35 (6.30-6.39) 


0.01 (<0.08) 


0.24 (0.09-0.44) 


61/42 


UGC 05101 


0.4-8.0 


1.0 (<6.5) 


1.9 (0.80-4.7) 


6.32 (6.27-6.44) 


O.Olf 


5.9 (2.7-630) 


15/8 


NGC 6240 


0.4-8.0 


0.4 (<1.9) 


1.4 (0.9-1.9) 


6.46 (6.39-6.57) 


0.20 (0.10-0.36) 


2.4 (1.7-3.2) 


217/129 



Note. — Fits to the 2.0-8.0 kcV bandpass arc based on the model given in Table 5 with the addition of a Gaussian component with the parameters 
shown here for the second (i.e., intrinsic) absorber, power-law slope and Gaussian line. Fits to the 3.0-8.0 kcV bandpass arc based on an absorbed 
power-law plus Gaussian line model. 
* Second observation 



Table 7. Chandra Fe-K Line Limits 



Table 8. Double Plasma + Power-law Fits to Mkn 273 Spectra 



Galaxy 


F Fo _ K (10~ 14 ergs s" 1 ) 


Method 


Arp 220 


<0.92 


Direct measurement 


IRAS 05189 


<5.3 


Spectral fit 


IRAS 17208 


<0.83 


Direct measurement 


IRAS 20551 


<0.88 


Direct measurement 


IRAS 23128 


<1.3 


Direct measurement 


Mkn 231 


<5.0 


Spectral fit 


Mkn 273 


6.0 (2.2-11) 


Spectral fit 


UGC 05101 


3.5 (1.6-380) 


Spectral fit 


NGC 6240 


30. (22.-40.) 


Spectral fit 



Parameter 



Nuclear Spectrum Global Spectrum 



Method used to determine the Fe-K flux, as described in the 
text. Briefly, for sources with insufficient counts for spectral fitting, 
"direct measurement" refers to determining an upper-limit to the 
Fe-K flux based on the background level at 6.4 keV in the rest 
frame. 
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Fig. 4. — Long-term light curves of ULIRG sample. 



05101 in the starburst-ULIRG category), bearing in 
mind that this empirical classification may not reflect 



N H i (10 20 cm" 2 ) 
kTi (keV) 
Zl (JZq) 
K-l* 

N H ,2 (10 22 cm" 2 ) 
kT 2 ' (keV) 
Z 2 (/Z Q ) 
K 2 * 

N H 3 (10 23 cm" 2 ) 
Line E (keV) 
Line a (kev) 
Line N+ 

r 

N (xlO 4 ) 
X 2 /dof 



0.5 (< 1.8) 
0.84 (0.68-0.97) 
l.Of 

0.75 (0.50-1.4) 
1.2 (0.7-1.9) 
1.7 (1.0-2.0)* 
l.Of 

12 (8.4-20) 
3.5 (3.0-4.3) 
6.35 (6.30-6.39) 
(< 0.08) 
7.2 (2.6-11.8) 
1.8f 

7.4 (5.9-9.6) 
41.2/41 



2.1 (<7.8) 
0.66 (0.59-0.71) 
0.15 (0.10-0.23) 
17 (1.1-31) 

1.3 (0.9-2.2) 

1.1 (0.87-1.4) 
0.9 (0.1-5.0)' 
22 (4.5-64) 

3.4 (2.7-4.2) 
6.37 (6.30-6.42) 
0.0 (< 0.10) 

7.2 (2.2-13) 
1.8f 

7.6 (5.7-10) 
118/103 



Note. — The model fitted to the spectra was: 

Abs{N H ^)[Mekal{kT 1 ,Z u K 1 ) + Abs{N H , 2 ) X Mekal(kT 2 , Z 2 , K 2 ) + 
Abs(N H [ 3 ) X (PL(r,N) + Gaussian(E,cr, N Lina ))], with Abs = 
absorption, Mckal— plasma, PL — power law and Gaussian — Gaussian 
components. 

*Plasma normalization in units of 10~ 9 ^ n f n ^f V 

density in cm -3 , n p = Hydrogen density in cm" 3 

distance to source in cm. 
^Gaussian line normalization in units of 10 — 6 photons cm -2 
* Parameter reached pre-set boundary during error search. 



= electron 
D = luminosity 



the true nature of the dominant energy source. For ex- 
ample, it is fairly clear that an AGN contributes sub- 
stantially to the hard X-ray emission from NGC 6240 
ijlwasawa fc Comastril 119981 IVignati et al.l 119991 ) even 



though it is an optical LINE R that appears starburst- 
dominated in the infrared _|Imanishi fc Dudley! 120001: 



ILutz. Veilleux. fc Genzel 1 1999(1 . We will see below that 
the AGN- and starburst-ULIRGs do differ significantly 
in their hard X-ray properties. 

6.1. Morphology 
In every case, significant hard X-ray emission is present 



in the nuclear region, 
clei are within 0.5-1.5 



In general, the hard X-ray nu- 
of radio and CO nuclear posi- 
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Table 9. Complex Spectral Fits to NGC 6240 Spectra 



Parameter 




Nuclear Spectrum 


Global Spectrum 


JMh,i (10 cm 


) 


0.4 (<5.5) 


0.5 (0.3-0.9) 


Nh,2 (10 21 cm" 


2 ) 


5.4 (<9.4) 


0.3 (<1.5) 


Kl 1 (KCV ) 




U.zo (U.l t-yj.oo ) 


c\ Qfi ( ' c\ o^* n qq\ 
U.oU {\j.lt-\j.6d) 


Z\ { Zq) 




2.6 (0.1-5.0)+ 


0.32 (0.11-1.1) 


* 




0.51 (0.04-1.9) 


4.1 (1.3-13.) 


iVfr,3 (io 22 cm" 


~ 2 ) 


1.3 (0.77-1.4) 


0.78 (0.69-0.90) 


fcT 2 ' (keV) 




0.67 (0.61-0.76) 


0.66 (0.62-0.72) 


Z2 {/Zq) 




9.0 (3.1-10)* 


0.76 (0.32-1.23) 






1.3 (0.88-3.9) 


21 (3.5-41.) 


Linel E (keV) 




2.31 (2.25-2.35) 


2.33 (2.31-2.35) 


Linel <r (kev) 




O.Olf 


O.Olf 


Linel N+ 




5.4 (2.0-9.0) 


12. (7.2-16) 


JVh,4 (10 22 cm- 


- 2 ) 


1.5 (0.5-2.3) 


1.1 (0.7-1.7) 


Gamma 




1.55 (0.98-2.02) 


1.62 (1.22-1.93) 


AT (xlO 4 ) 




1.8 (0.9-3.5) 


2.7 (1.3-4.2) 


Line2 E (keV) 




6.35 (6.33-6.38) 


6.34 (6.31-6.37) 


Line2 <r (kev) 




0.00 (<0.06) 


0.02 (<0.08) 


Line2 N+ 




16 (11-21) 


17 (11-24) 


Line3 E (keV) 




6.57 (6.50-6.69) 


6.62 (6.54-6.67) 


Line3 a (kev) 




0.01 (< 0.20) 


0.01 (<0.13) 


Line3 N+ 




4.9 (0.8-11) 


7.4 (2.2-13) 


Nufi (10 24 cm- 


- 2 ) 


1.3 (0.5-3.0) 


1.3 (0.7-2.3) 


N re fl 




2.5 (0.5-5.7) 


4.6 (1.4-13) 


X 2 /do/ 




147/117 


331/225 



Note. — The model fitted to the spectra was: Abs(Nn,\) X 
[Abs(N H . 2 ) x Mekal{kT x , Z t , K 1 )+Abs(N H , 2 )x [Mekal(kT 2 , Z 2 , K 2 ) + 
GaussianiE!^!, N LincA ) + Abs(N H:3 ) X [PL(T, N) + 
Gaussian{E 2 , <r 2 , Nhine, 2 ) + Gaussian(E 3 , CT3, N Line,s)\ + 
Abs(N H i ) X Refl(N rl! fi)], with Abs = absorption, Mckal=plasma, PL 
— power law, Gaussian — Gaussian components and Ref! — Compton 
reflection components. 

*Plasma normalization in units of 10~ 10 ^ "e"H dy ' n — electron 
density in cm -3 , n p = Hydrogen density in cm -3 , D = luminosity 
distance to source in cm. 

t Gaussian line normalization in units of 10 — 6 photons cm -2 s — 1 . 

* Parameter reached pre-set boundary during error search. 



tions (Brvant & Scovillc 1999: Condon ct al. 1991, 


1996 


Planesas. Mirabel. & Sandersl 11991 


Scovillc et al. 


1991 


Thean et alJ l200Ct lYun & Scovilld 


1995). However, as 



noted above, in three cases (the starburst-ULIRGs Arp 
220, and IRAS 17208-0014) and in NGC 6240, this com- 
ponent is resolved on arcsecond (~ kpc) scales (see also 
iKomossa et aLlfeOOffi . In the other six cases where the 
nuclear source is unresolved, the upper limit to the size 
(~ 0.5 kpc FWHM) is consistent with either an AGN or 
compact starburst with a size similar to the nuclear radio 
continuum and CO emission (see Table Upland references 
therein) . 

In addition to the nuclear emission, a somewhat sur- 
prising result is that often extended hard X-ray emis- 
sion is detected beyond the nuclear region, on scales of- 
order ten kpc. In five cases (the starburst ULIRGs IRAS 
17208-0014, IRAS 20551-4250, IRAS 23128-5919, UGC 
5101) and in NGC 6240, this large-scale emission com- 
prises a significant fraction (~ 10 to 50%) of the total. 
Thus, in the five starburst-ULIRGs, the nuclear source 
is spatially resolved (two cases) and/or a significant frac- 
tion of the total hard X-ray emission arises well outside 
the nucleus (four cases) . The nuclear source is unresolved 
and dominates the total hard X-ray emission in the three 
AGN-ULIRGs. 

By way of comparison, Chandra observations of nearby 



starburst galaxies such as M82 and NGC 253 have 
shown that a large fraction of the hard X-ray emis- 
sion in these galaxies is due to point-sources, most 
likely high-mass X-ray binaries, concentrated within the 
central ~ kpc of the galaxy. Several of the ULIRGs 
have extra-nuclear sources that would qualify as ULXs 
if they are not interlopers. ULXs are often observed 
in starburst galaxies and, more notably for compari- 
son with ULIRGS , in m ergers such as the Antennae 
l|Zezas fc Fabbianol 12002]) . In IRAS 20551-4250 - 50% 
of the extra-nuclear hard X-ray flux may be due to a sin- 
gle ULX. The remainder of the hard X-ray emission in 
starbursts is due to an u nresolved population of lower- 
lumin osity X-ray binaries ijGrimm. Gilfanov. &: Sunvaev 
1 20031), very hot gas (> few keV, see iCappi et al. 
1999; Griffith s~et~atl l200fl). or inverse-Compton scatter- 
ing of IR photons off of relativistic electrons generated 
by supernova shocks llMoran. Lehnert. fc Hclfand Tl999t 
iPersic fc Rephaelill2002D . A strong point source at the 
dynamical center of the galaxy (i.e. a plausible AGN can- 
didate) is typically not detected (although cf., Weaver et 
al. 2002). 

The hard X-ray morphology in these ULIRGs is gen- 
erally consistent with these types of starbursts placed at 
the same distance as the ULIRG galaxies (even for the 
AGN-ULIRGs) . Table El also lists the semimajor axis 
size of nuclear CO observed in several ULIRGs along with 
the hard X-ray semimajor axis derived from Gaussian x 
PSF fits to the hard X-ray surface brightness. These val- 
ues are clearly consistent with the X-ray emission origi- 
nating from within the nuclear molecular cloud/disk, or 
from a region of similar size (in the cases of Arp 220 and 
IRAS 17208-0014). 

6.2. Spectroscopy 

As stated above, most of the galaxies in this sam- 
ple require at least two model components to fit their 
X-ray spectra. The soft X-ray spectra of these galax- 
ies show that hot gas with kT ~ 0.7 keV is present 
within the nuclei of these galaxies (with the possible 
exception of IRAS 05189-2524 where the AGN may be 
completely swamping any thermal flux). In this regard, 
they are very simila r to typical starburst galaxies ob- 
served with Chandra iLira , Ward. Zezas Murrav^Q Q 2[ 
Martin. Kobulnickv. fc HecbnsLu|l27l02TlSt 
2002). 

The hard X-rays we observe may in principle be due to 
either an AGN or starburst. Our analysis shows that the 
hard X-ray emission detected in the ULIRGs is generally 
not absorbed by large column densities (i.e., in excess 
of 10 23 cm" 2 ). Mkn 273 is the only exception (see also 
iXia et aUl200l . If the hard X -rays are due to an AGN, 
then there are two possibilities: the AGN is truly unab- 
sorbed (and hence has a low intrinsic luminosity) or the 
AGN is "Compton-thick" (i.e., absorbed by a column in 
excess of 10 24 cm -2 ) and the observed hard X-rays are 
scattered into the line of sight. In the first case we would 
expect the hard X-ray slopes to be similar to those ob- 
served in Seyfert 1 a nd Compton-thin Seyfert 2 galaxies , 
namely T = 1.5-2.0 Jbahlem. Weaver & Heckmanll998t 
iMushotzkv. Done, fc Pounds! 11993: iPtak et alJ 11999(1 
IRAS 05189-2524 is consistent with this perspective, once 
we allow for the possible presence of Fe-K emission with 
an EW of ~ 0.1 keV (see TableEJl- On the other hand the 
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Table 10. Fits to Global Spectra 



Parameter 


Arp 220 


105189 


117208 


120551 


123128 


Mkn 231 


UGC 05101 


IV J l^lu Cm j 


n 9n +0.94 | 
U - ZU -0.00 


, oq+0.48 
— 1 25 


c- 99 +4.02 
°- zz -4.72 


n on+ 2 ^ 2 t 


9 QS + 2, 17 

2.y»_ 1 n 


n n^+o.79 j 


9 9Q+3-12 
z - za -2.14 




n on+0.06 

u - oo -0.03 


30+ ' 31 * 
u - ou -o.oo 


35+ ' 96 


32+ 012 

u '° -0.02 


50+° 19 


53+ 15 

u - do -0 .23 


70+ - 34 

u - lu -0.40 


^i (/Za) 


l.Of 


„ „+n t 

0-0-ao * 


^ -. +4 9 

O-i-o.i 


l.Of 


^ -1 -1-0 1 

O-l-o.o 


„ „+0 3 
- 2 -0i 


^ -. +0 l 

0-l-o.o 


V * 

Ki 


1.20 


141.40 












AT ( i n22 „»~, — 2\ 

^if,2 (10 cm J 


n 7+0-2 

0-7-0.5 


in c+89.4 
10.6_ g g 














0.8_ 2 


9 n+0-0 1 












^2 (/Z ) 


l.Of 


6+ 4 ' 4 
u '°-0.5 












^2* 


8.84 


563.70 












^H,3 (10 22 cm" 2 ) 


1 0+ 29 

!- U -1.0 






0.7181 


e 1+6.8 
-!-3.9 


0+ 2 1 * 


7+ 4 - 6 


r 


1 14+0-93 
!•! -0.87 


0.52«;f 3 




0.85± 2 ;° 4 


2 5Q+ 1 - 95 


u -°°-0.41 


l 06+ 1 - 79 
1 - UD -1.01 


^3 (X10 5 ) 


1.38 


11.32 


1.53 


0.51 


21.59 


2.85 


0.80 


Line E (keV) 














6 Q0+ - 00 
o.au_ 29 


Line a (keV) 














o 7 +0.13 
u - 31 -0.37 


Line N+ (xlO 6 ) 














a co+68.60 
4.D0_2 gl 


X 2 /dof 


48.9/64 


94.4/75 


18.7/16 


21.0/27 


30.7/30 


83.1/89 


32.5/16 



* Plasma normalization in units of 10 9 ", " "i5> — 1 «e = electron 

density in cm -3 , n p = Hydrogen density in cm -3 , D = luminosity 
distance to source in cm. 

^Gaussian line normalization in units of 10 — 6 photons cm -2 s — 1 . 

t Parameter reached pre-set boundary during error search. 



Table 11. Fluxes and Luminosities 



Galaxy 




Global 






Nuclear 


0.5-2.0 keV 


2-10 keV 


-?b. 5-2.0 keV 


F2-IO keV 


-^2-10 keV 


F0.5- 


-2.0 keV -F2-10 keV 


F N 

Fg 


Cn 
C G 


F N 

f g 


Cg 


Arp 220 


0.6 


1.3 


0.10 


0.16 


1.2 


0.27 


0.25 


0.85 


0.65 


IRAS 05189-2524* 


0.73 


37. 


15. 


0.66 


35. 


0.90 


0.87 


0.96 


0.95 


IRAS 17208-0014 


0.21 


0.62 


0.25 


0.11 


0.42 


0.52 


0.43 


0.68 


0.77 


IRAS 20551-4250 


0.49 


0.84 


0.35 


0.18 


0.39 


0.37 


0.37 


0.46 


0.42 


IRAS 23128-5919 


0.43 


1.5 


0.66 


0.13 


1.3 


0.31 


0.25 


0.92 


0.67 


Mkn 231 


1.0 


7.8 


3.1 


0.53 


8.1 


0.53 


0.47 


1.04 


0.85 


Mkn 273 


1.1 


8.3 


2.6 


0.34 


7.8 


0.30 


0.26 


0.95 


0.88 


UGC 05101 


0.20 


1.6 


0.55 


0.12 


0.81 


0.60 


0.52 


0.51 


0.76 


NGC 6240 


7.2 


26. 


3.4 


1.8 


17. 


0.26 


0.23 


0.66 


0.63 



NOTE. — Fluxes arc in units of 10 13 erg cm 2 s 1 , luminosities arc in units of 10 42 erg s , derived from plasma plus power-law fits. The 
"nuclear" 2-10 keV flux exceeds the "global" flux in Mkn 231 due to statistical uncertainty in the power-law slope; the global 2-10 keV flux of Mkn 
231 is 8.9 X 10 — 13 erg cm -2 s _1 whcn the power-law slope is fixed at best-fit value of the nuclear spectrum fit (T — 0.21). is the ratio of nuclear 

and global fluxes while S^- is the ratio of nuclear and global count rates. 



Table 12. ASCA and BeppoSAX Fluxes and Fe-K Line Limits 



Galaxy 


Satellite 


Date 


Fo. 5-2.0 keV 


^2-10 keV 


Fe-K Line EW (keV) 


References 


Arp 220 


ASCA 


03/27/1994 


0.6 


1.8 




1 


Arp 220 


BeppoSAX 


08/04/2000 




2.2 




7 


Arp 220 


BeppoSAX 


08/27/2000 




1.7 




7 


Mkn 231 


ASCA 


12/05/1994 


1.2 


6.0 


< 0.9 


1 


Mkn 273 


ASCA 


10/27/1994 


1.6 


5.6 


0.52 (0.26-1.0) 


1 


Mkn 273 


BeppoSAX 


01/25/1999 




3.5 


1.2 (0.2-3.2) 


2 


IRAS 05189-2524 


ASCA 


02/15/1995 




48 


0.11 (0.03-0.18) 


4 


IRAS 05189-2524 


BeppoSAX 


10/03/1999 




36 


0.14 (0.05-0.34) 


1 


IRAS 17208-0014 


ASCA 


03/06/1996 




3.0 




2 


IRAS 20551-4250 


ASCA 


10/19/1995 




3.1 




3 


IRAS 23128-5919 


ASCA 


02/11/1995 




3.4 




3 


UGC 05101 


ASCA 


11/17/1996 


1.7 


6.2 




7 


NGC 6240 


ASCA 


03/27/1994 


6.4 


19. 


2.2 


5 


NGC 6240 


BeppoSAX 


08/15/1998 


7.4 


32. 


1.6 


6 



References. — 1 llwasawal I1999D: 2. IRisaliti. Gilli. Maiolino. &: SalvatH 120001) : 3. IMisaki et"aT1 I1999D : 4. ISevergnini et al.l 120011) : 5. 

Iwasawa &c Comastri I 1998); 6. Vignati ct al .1 11999ft : 7.' This work 



Note. — Fluxes in units of 10 13 erg cm 2 s 1 . 
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Table 13. ULIRG Parameters and Classification 



Galaxy 


Lfir 




^2-10 kcV 

Ffir 


Classification 


log N Ht co* 


Rco ] 


-^2-8 kcV* 




(1O 12 L ) 


(io- 4 ) 


Optical 


IK 


(cm -2 ) 


(") 


(") 


Arp 220 


0.9 


0.08 


0.27 


LINER 1 


HII 2 


24.7 


0.9 


0.9 (0.8-1.1) 


IRAS 05189 


0.7 


0.25 


56 


S2 3 


AGN 2 , SI 4 






< 0.28 


IRAS 17208 


1.6 


0.05 


0.44 


HII 1 


PDR 5 


24.1 


1.25 


1.0 (0.7-1.2) 


IRAS 20551 


0.6 


0.15 


1.5 


HII 6 


HII 5 






< 0.70 


IRAS 23128 


0.6 


0.15 


4.8 


HII 6 


HII 5 






< 0.47 


Mkn 231 


1.6 


0.27 


8.6 


SI 3 


AGN 2 


24.6 


0.5 


< 0.25 


Mkn 273 


0.8 


0.11 


5.8 


S2 1 


AGN 2 


>24.6 


< 0.6 


< 0.47 


UGC 05101 


0.6 


0.09 


1.9 


LINER 7 


AGN 8 , SB 9 






< 0.44 


NGC 6240 


0.4 


0.15 


17 


LINER 7 


SB 9 


24.6 


0.39 


1.1 (1.0-1.2) 



Note. — Lfir — 2.58Leo/im 4- 

Column density (i n cm" 2 ) inferred from t he CO mass surface 
brightnesses given in Bryant & Scovillc (1999). 

^Deconvolved semima jor axis of CO emission given in 
IBrvant fc Scoville! tl999l) . 

' Deconvolved semimajor axis of 2-10 keV emission, based on el- 
li ptical Gaussian fits (includ ing PSF model) discussed in the text. 

1Kirm_\^meux 1 _fcS^^ 

: Imanishi & Dudlcv 12000) 

: Vcilleux, Kim, k. Sanders (1999) 

4 Vcillcux, Sanders. & Kim (1999) 

: Laurent eTaTI 1^200(11 

{ Kewlev^e^er^o^ita^^Lurnsde^ Jjffli^ 

l^nleuj^Km^arid^r^Mazza^lla^^orfen 119951) 

6 Imanishi. Dudley & Malonev (2001) 

E Lutz. Vcillcux. & Gcnzcl ( 1999) 
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hard X-ray spectrum of Mkn 231 appears to be genuinely 
flat (r < 1). 

In the Compton-thick scenario for an AGN origin 
to the hard X-rays, the hard X-ray would be due to 
some combination of reflection from optically-thick neu- 
tral material, scattering from optically-thin (and likely 
highly-ionized) material and leakage of X-rays through 
patches in the obscuring material. These effects would 
tend to flatten the observed X-ray spectrum and pro- 
duce high-EW (> 1 keV) Fe-K lines, as observed in 
samples of Sevfert 2 llTurner et al.lll997l) and composite 
starburst-Seyfert l)Levenson. Weaver, fc Heckmanl 12001^ 
galaxies (note that scattering from highly-ionized ma- 
terial alone would not flatten the observed X-ray spec- 
trum). The fact that the hard X-ray emission tends 
to be coincident with CO gas with high implied col- 
umn densities (see Table EH supports this picture (al- 
though a caveat is discussed below), and this is clearly 
the case in NGC 6240. The Fe-K line EW in Mkn 273 
is consistent with that expected to be produced by the 
transmission of X-rays throug h material with colum n 
densities of order 10 cm" " 2 lILeahv fc CreightonllTOQ^ . 
as is observed directly. However, there are insufficient 
counts in most of these spectra to allow meaningful con- 
straints to be placed on complex models (i.e., involv- 
ing reflection and "leaky" absorbers). Partial-covering 
fits to the XMM-Newton spectra of IRAS 05189-2524 
and UGC 05101 and the Chandra spectrum of Mkn 
231 (|Oallagher et alJ [2002) are in fact consistent with a 
highly-obscured hard X-ray source that is leaky or con- 
tains one or more scattering or reflection regions. Similar 
results based on XMM-Newton o bservations of ULIRGs 
are given in IBraito et alJ l|2002() . Also note that if a 
highly-ionized scattering medium is acting as a mirror 
for the hard X-rays, then the Fe-K lines expected from 
such a mirror would be due to Fe XXV and Fe XXVI 
(Ka energies o f 6.7 and 6.9 keV) with EWs on the order 
of 0.5-1.0 keV (Bi anchi fc Matdl200^ . We suggest that 
in these AGN ULIRGs the amounts of reflected and/or 
leaked flux and scattered flux are comparable to each 
other, effectively reducing the equivalent widths expected 
from neutral and ionized Fe-K to values within the range 
of our upper- limits. In the case of IRAS 05189-2524, the 
very tight constraint on neutral Fe-K suggests that the 
flux is mostly scattered and not reflected (or again IRAS 
05189-2524 is simply Compton-thin). 

Note that the long-term variability observed in the 
AGN ULIRGs implies that both the scattering and re- 
flecting material must be located within 1-2 pc of the 
central X-ray sources, or distances consistent with the 
expected size the putative tori. However, this argues 
against the arcsecond-scale CO gas being a source of re- 
flection that dominates the hard X-ray flux in the AGN 
ULIRGs. This is consistent with the relatively high 
£2-10 kcv/Lm ratio of AGN ULIRGS, i.e., the molecular 
gas is most likely only directly obscuring the AGN with 
columns of order 10 23 cm~ 2 or less. In general the CO 
gas appears to be distributed in 100-500 pc scale disks, 
and the disks are either viewed face-on and/or (partic- 
ularly in the case of double n uclei) the nuclei are not 
located centrally in the disks iJBrvant fe Scovillel 119991 
IDownes fe Solomon! Il998|) . Both of these effects would 
result in lower column densities being observed toward 
the nuclei than the peak values derived from the CO 



data. 

Turning to the starburst-ULIRGs, the hard X-ray 
emission is too weak for strong constraints to be placed 
on its spectral form (it is consistent with either a weak 
AGN viewed through modest columns, or a typical star- 
burst). Moreover, the upper limits on the Fe Ka EW are 
not restrictive. The extended hard X-ray emission re- 
solved by Chandra in NGC 6240 suggests that starburst 
processes such as these may be contributing significantly 
(i.e., the two nuclei contribute <50% of the 2-10 X-ray 
flux o f NGC 6240 although, as shown in Komoss a et alJ 
(2003), the nuclei dominate in the 5-8 keV band). These 
possibilities will be discussed in more detail in a subse- 
quent paper. 

6.3. The X-ray Emission of ULIRGs in Context 
6.3.1. The Far- Infrared 

We have argued above that the signature of an AGN is 
possibly present in the hard X-ray emission of at least 
three of the eight ULIRGs in our sample (the same 
three that show evidence for an AGN in their optical 
and IR spectra), and that Fe-K emission shows that an 
energetically-important AGN is present in two of the 
starburst ULIRGs. One way to assess possible energetic 
contribution of AGN is to compare the ratio of the hard 
X-ray to far-infrared luminosities in the ULIRGs to val- 
ues in typic al AGN and starbursts. 

Following iLevenson. Weaver, fe Heckmanl 1)20011 here- 
after LWH), we first plot the ratio of hard X-ray 
to FIR luminosity as a function of the i^Sjum/ 'Feoum 
flux ratio (Figure [SJ. The latter (a measure of 
luminosity-weighted mean dust temperature) is a use- 
ful empirical diagnostic to help assess the relative 
energetic importa nce of AGN and starburst activ- 
ity in a galaxy ijde Griip. Milev. fe Lubl Il987j) . We 
compare the ULIR Gs to a sample of starburs t s ob- 
served by ASCA llDahlem Weaver fc Heckmanl 119981: 
[ Delia Ceca et al.ll996tll999tlMoran Lehnert. fe Helfandl 
Il999t iPtak et al.lll999t iTurner et alJll997|) . a sample of 
composite starburst /Seyfert 2 galaxies, and the "pure" 
Seyfert 2 and Seyfert 1 galaxies samples from LWH. 
As shown by LWH, the pure Seyfert 2's are generally 
Compton-thin and the composite starburst /Seyfert 2's 
are generally Compton-thick. Since these other samples 
are based on ASCA data, the ULIRG 2-10 keV luminosi- 
ties were taken from the global spectral fits, although we 
also plot the points derived from the nuclear fluxes for 
comparison. The IRAS infrared fluxes have been taken 
from NED, and starburst galaxies for which the global 
fluxes fluxes differ by more than 10% from the poin t- 
source values have been taken from lSoifer et alJ {l989). 

Several results are clear from this comparison. First, 
both the mid/far-IR color and the relative strength of the 
hard X-ray continuum emission in the ULIRGs is gen- 
erally similar to the values for starburst galaxies. The 
latter ratio is typically about two orders-of-magnitude 
below the values for most type 1 Seyferts and pure 
(Compton-thin) type 2 Seyfert galaxies. IRAS 05189- 
2524 is the only ULIRG lying outside the starburst 
regime and overlapping the pure type 2 Seyfert regime. 
Note that we are using observed hard X-ray luminosi- 
ties in these plots, and correcting for the ~ 10 23 cm~ 2 
column density in Mkn 273 would increase its 2-10 keV 
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Fig. 5. — 1-2-10 kcv/LiR vs. F2 5fl /F 60M . The plot key is: 
ULIRGS=red circles (values shown with filled circles are based 
on X-ray fluxes derived from the nuclear spectrum), star bursts = 
green filled circles, Seyfert 2s = blue squares, composites = cyan 
filled squares, Seyfert Is = purple diamonds. 



luminosity by a factor of ~ 3 and place it in the HX/FIR 
> 10~ 3 Compton-thin regime. There is some overlap be- 
tween the ULIRGs and the (Compton-thick) composite 
starburst/Seyfert 2 galaxies, but the ULIRGs as-a-class 
are even more extreme. Thus, in a purely empirical ener- 
getic sense there is a need to invoke a significant contri- 
bution of an AGN to the hard X-ray emission in only one 
of the eight ULIRGs. Second, the three AGN-ULIRGs 
(and NGC 6240) have significantly higher ratios of hard 
X-ray to far-IR luminosity than the starburst-ULIRGs. 
This suggests that in fact an AGN is contributing signif- 
icantly to the hard X-ray emission in these former cases. 

The above conclusions are reinforced by Figure HJI 
where we have plotted the ratio of Fe Ka and far-IR 
luminosities vs. the ratio of hard X-ray and far-IR lumi- 
nosities for the above samples of ULIRGs and AGN. In 
particular, the upper limits to the ratio of the Fe Ka and 
far-IR fluxes in the ULIRGs lie below even the very low 
ratios observed in the Compton-thick starburst/Seyfert 2 
composites, and are up to two orders-of-magnitude lower 
than in the pure type 2 Seyferts. 

6.3.2. The [OIIIJX5007 Line 

The [OIII]A5007 emission- line is the strongest optical 
line produced in the kpc-scale Narrow Line Region, and is 
often used as a rough indicator of the the true luminosity 
of the AGN in both type 1 and type 2 Seyferts. On this 
basis, Bas sani et all l)1999|) have proposed using the ratio 
of the hard X-ray and [OIII]A5007 fluxes as an indicator 
of X-ray absorption that is valid even in the Compton- 
thick regime. They also show that there is a strong in- 
verse correlation between the hard X-ray/ [OIII] flux ra- 
tio and the Fe Ka equivalent width, with Compton-thick 
type 2 Seyfert nuclei at one extreme and type 1 Seyferts 
at the other. 

In Figure [7| we plot the luminosities of the [OIIIJA5007 
lines vs. those of the Fe Ka lines for Seyferts and 
ULIRGs. We have also corrected the [OIII] for dust 
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Fig. 6. — L 2 -io kev/^/fl vs - L Fe - K /L IR . The upper-limits 
are based on fluxes listed in Table 6. 



extinction using th e prescription (based on the Balmer 
decrement) given in Dahari & D e Roberts! (|L088|) . In the 
case of the ULIRGs we used the Ha/H/3 flux ratio listed 
in TablelTHwhile we applied mean corrections (also based 
on the statistical analysis in iDahari fe De Robertsll988|) 
to the Seyfert 1 and Seyfert 2 samples (which increased 
the [OIII] fluxes by factors of 4 and 10, respectively). 
In the plot of the extinction-corrected [OIII] luminosity, 
the Seyferts exhibit a correlation with roughly unit slope 
and a scatter of ~ ±0.5 dex. This suggests that both 
the Fe Ka line and the [OIII]A5007 line can be used as 
very rough indicators of the luminosity of the hidden 
AGN, even in Compton-thick type 2 Seyferts. Since we 
have Fe Ka detections for only two of the eight ULIRGs, 
we can only say that the data are consistent with the 
ULIRGs following the same trend as the Seyferts. More 
to the point, the relative weakness of the Ka line in the 
ULIRGs (Figure HJJl suggests that energetically dominant 
AGN are not present in most ULIRGs (even though AGN 
may contribute significantly to the hard X-ray emission 
in some cases). 

To further assess the possible contribution of an AGN 
to the hard X-ray emission in ULIRGs, we have plotted 
the ratio of the hard X-ray and [OIII]A5007 flux (see Ta- 
ble ll4f> vs. the ratio of the hard X-ray and far-IR flux for 
our samples of ULIRGs and Seyferts (Figure [SJ. Prior 
to [OIII] extinction correction the ULIRGs occupy a dis- 
tinct part of parameter space in this figure. On the one 
hand, they have ratios of hard X-ray to [OIII] flux that 
are most similar to type 1 and pure (Compton-thin) type 
2 Seyfert galaxies, and are significantly larger than in the 
Co mpton-thick starburst /Seyfert 2 composites. Follow- 
ing Bassani ct al. (1999), this implies only modest X- 
ray absorption in the ULIRGs. On the other hand (as 
discussed above) the ULIRGs have extremely small ra- 
tios of hard X-ray to far-infrared flux. After extinction- 
correcting the [OIII] fluxes, the ULIRGs better overlap 
the locus of the Compton-thick Seyfert 2 and composite 
galaxies, albeit with much more dispersion. The larger 
dispersion in the ULIRG sample is due at least in part to 
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Fig. 7. — Fe-K luminosity plotted as a function of [OIII] line luminosity, (left) Plot generated using observed [OIII] luminosities; (right) 
plot generated using extinction-corrected [OIII] luminosities (see text). The plot key is: ULIRGS=red circles, Seyfert 2s = blue squares, 
composites = cyan rilled squares, Seyfert Is = purple diamonds. 



the large (and hence uncertain) corrections to the [OIII] 
fluxes. Taking the extinction-corrected [OIII] fluxes at 
face- value, the AGN-ULIRGs Mkn 231 and IRAS 05189- 
2524 and the starburst-ULIRG Arp 220 have hard X-ray 
to [OIII] flux ratios similar to Compton-thin Seyferts, 
the AGN-ULIRG Mkn 273, the starburst-ULIRG UGC 
05101, and NGC 6240 have ratios similar to Compton- 
thick Seyferts, and the starburst-ULIRGS IRAS 17208- 
0014, IRAS 20551-4250, and IRAS 23128-5919 are inter- 
mediate. 

7. SUMMARY 

We have presented the initial results from a Chandra 
survey of a complete sample of the 8 nearest (z < 0.04) 
ultraluminous IR galaxies (ULIRGs) plus NGC 6240, us- 
ing the hard X-rays (2-8 kev) to search for the possible 
presence of an obscured AGN, and to assess its contribu- 
tion to the bolometric luminosity. In six cases (includ- 
ing NGC 6240), the extant optical and infrared spectra 
suggest that a starburst dominates the energetics (the 
" starburst-ULIRGs" ) , while an energetically significant 
AGN is present in the other three (the "AGN-ULIRGs"). 
We find that the hard X-ray properties of these two sub- 
types differ as well. 

A hard X-ray source is detected in the nuclear region 
in every case. The nuclear source is spatially-resolved in 
two of the starburst-ULIRGs Arp 220 and IRAS 17208- 
0014 (FWHM ~ 1 kpc and 4 kpc respectively) and also in 
NGC 6240, and is unresolved (FWHM < 0.5 kpc) in the 
others. The upper limits to the spatial extent in the six 
unresolved cases are consistent with either an AGN or a 
compact starburst. Hard X-ray emission on larger (galac- 
tic) spatial scales is significant in six cases (all starburst- 
ULIRGs), comprising 10 to 50% of the total flux. We 
have shown that five starburst-ULIRGs have hard X-ray 
luminosities about an order-of-magnitude smaller than 
the three AGN-ULIRGs, with the hard X-ray luminosity 
of NGC 6240 being comparable to the AGN-ULIRGs. 



Our analysis of the hard X-ray spectra provides no di- 
rect evidence for absorbing columns in excess of ~ 10 23 
cm~ 2 except for in Mkn 273, although the data are rel- 
atively poor for the five starburst-ULIRGs. The Fe Ka 
line is convincingly detected in Mkn 273 and NGC 6240 
and marginally detected in UGC 05101. 

The ratio of the hard X-ray to far-IR flux in the 
ULIRGs (HX/FIR) is about three (two) orders of mag- 
nitude smaller than in type 1 (Compton-thin type 2) 
Seyfert galaxies. Only the AGN-ULIRG IRAS05189- 
2524 has a value for HX/FIR that is significantly higher 
than that in typical starburst galaxies, and only IRAS 
05189-2524 and Mkn 231 are in the Compton-thin regime 
of a HX/FIR vs. HX/[OIII] plot. The three X-ray- 
brightest ULIRGS (the AGN-ULIRGs) and NGC 6240 
have intermediate values of HX/FIR similar to the 
Seyfert 2/ starburst composite systems studied by LWH. 
Likewise, the flux ratio of the Fe Ka line and far-IR 
continuum is usually at least two orders-of-magnitude 
smaller in the ULIRGs than in typical Seyfert galaxies 
(but with some overlap with the Seyfert2/starburst com- 
posites) . We show that type 1 and type 2 Seyferts follow 
a correlation between the [OIII]A5007 and Fe Ka lumi- 
nosities. The weakness of the [OIII] and Fe Ka emission 
in the ULIRGs (relative to the far-IR) suggests a corre- 
spondingly weak AGN, although a highly-absorbed AGN 
(possibly) in conjunction with a complex scattering ge- 
ometry cannot be ruled out. 

Existing optical and IR spectroscopy show that AGN 
are definitely present in three of the ULIRGs, and we con- 
clude that these AGN make a significant contribution to 
the observed hard X-ray emission in these cases. There is 
no compelling reason to invoke the presence of an AGN 
in the other five cases except for the strong Fe-K emission 
in UGC 05101 and NGC 6240 (where the AGN contri- 
bution to hard X-rays is also well established from the 
BeppoSAX PDS data; Vignati et al. 1999). Even in the 
three AGN-ULIRGs, the contribution of the AGN to the 
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Table 14. ULIRG [OIII] Data 
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which dominates the hard X-ray emission. 
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Fig. 8. — 2-10 keV/FIR luminosity plotted as a function of 2-10 keV/[OIII] luminosity, (left) Plot generated using observed [OIII] 
luminosities; (right) plot generated using extinction-corrected [OIII] lumi nosities ( see te xt). The box delineated in the upper-right corner 
shows the region where Compton-thin Seyfert 2s are found in[Pancssa & Bassanil Pool) . The plot key is: ULIRGS=red circles, Seyfert 2s 
= blue squares, composites = cyan filled squares, Seyfert Is = purple diamonds. 



bolomctric (IR) luminosity is quite uncertain. To make 
a major contribution, an AGN must be buried behind 
highly Compton-thick material with a very small trans- 
mitted or reflected fraction compared to typical Seyfert 
galaxies. This could be possible (given the very large 
column densities of molecular gas observed in the nuclei 
of these galaxies), but our new data provide no direct 
evidence for such high absorbing columns. 

Of course, absence of evidence is not evidence of ab- 
sence: it is exceedingly difficult to robustly prove that 
powerful AGN can not be present in these galaxies. How- 
ever, we conclude that our new data provide no evidence 
that powerful "buried quasars" dominate the overall en- 



ergetics of most ULIRGs, particularly those with a star- 
burst optical or IR classification. 
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